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Morphological variation of A/burnus mossulensis (Heckel, 1843) populations in Iran were 
investigated by collection of 705 specimens from 27 rivers of five basins, including 
Bushehr, Fars, Karkheh, Karun and tributaries of Tigris (Diyala) basin in Iran in 2010 by 
a seine net. The specimens were fixed in 10% formalin for further investigation after 
anesthetizing in 1% clove oil solution. Twenty-two morphometric and 11 meristic 
characters were examined. Morphometric characters in adjusted form and meristic 
characters in classified form were used for population comparison. Analysis of 
morphometric and meristic characters by ANOVA showed significant differences 
(p<0.05) in all basins for all characters except the number of dorsal and anal fin spines. 
Classification of meristic characters showed most specimens of all basins have eight soft 
rays in dorsal fin, 12 soft rays in anal fin, 19 soft rays in caudal fin, 15 soft rays in pectoral 
fin and nine soft rays in pelvic fin. Discriminant function analysis showed that 
populations from Fars and Karun Basins were different from each other and from other 
populations, on the other hand, populations from Diyala, Karkheh and Bushehr basins 
overlapped. This result may indicate that there are similar conditions in Diyala, Karkheh 
and Bushehr basins, resulting in more similarity amongst the populations of these basins. 
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INTRODUCTION 

Investigation of aquatic ecosystems, especially fishes are important from view point of evolution, 
ecology, behavior, conservation, stock identification, stock assessment of fishes and water resource 
management (Anvarifar e¢ a/, 2011). The study of morphological characters with the purpose of 
identification or characterizing fish stock units for some time has been a valuable interest in 
ichthyology (Tudela, 1999). The study of morphometric and meristic characters are forceful tools for 
measuring and differentiation of the same species (Naeem & Salam, 2005). A sufficient amount of 
separation may result in notable morphological, meristic and genetic differentiation among stocks 
within a species, which may be recognizable as a basis for the management of different stocks (Mian 
et al, 2014). Meristic characters are most commonly used for differentiation of species and 
populations (Sedaghat ef a/, 2012). Morphometric measurements and its statistical relationships of 
fishes are essential for both fishery biology and taxonomic studies (Gharaei, 2012; Poria et a/., 2013; 
Sharma e¢ a/., 2014; Zamani-Faradonbe ef a/, 2015b). Investigation of morphometric and meristic 
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characters remains the simplest and most direct way among methods of species identification. It is 
well-known that the examination of phenotypic differences in morphometric characters or meristic 
counts is the technique most commonly used to define stocks of fish (Samaradivakara et a/., 2015; 
Keivany & Arab, 2017; Mouludi-Saleh & Keivany, 2018). Despite the introduction of techniques 
which directly examines biochemical or molecular genetic variation, these conventional methods 
continue to have an important role in stock identification even to date (Swain & Foote, 1999; Jalili e¢ 
al, 2015; Keivany ef a/, 2015). 

The cyprinid genus, A/burnus Rafinesque, 1820 comprises about 39 recognized species 
distributed in Europe and West Asia (Falahatkar ef a/, 2015). This genus has eight confirmed species 
in Iranian inland waters (Keivany e¢ a/., 2016d), including A. chalcoides (Gildenstadt, 1772), A. filippii 
Kessler, 1877 and A. hohenackeri Kessler, 1877 in the south Caspian Sea basin, A. atropatenae Berg, 
1925 in the Urmia Lake basin, A. mossulensis Heckel 1843 in the Tigris River, Esfahan, Kor River, 
Lake Maharlu, Bushehr and Hormuz basins, A. zagrosensis Coad, 2009 and A. caeruleus Heckel, 1843 
in the Tigris River basin and A. amirkabiri, Mousavi-Sabet ef a/., 2015 in the Namak Lake basin 
(Mousavi-Sabet e¢ a/., 2014, 2015). In addition, A. doriae de Filippi, 1865 and A. maculatus Keyserling, 
1861 have uncertain provenance and validity from Iran. Determination of fish stocks is too 
important for rational and effective fisheries management, because each stock needs separate 
management to aim of optimal harvest (Salini e¢ a/, 2004; Erguden & Turan, 2005). Despite some 
works on the biology of this species (Ketvany e¢ a/, 2016a, b, c; 2017a, b), there is no thorough 
examination of the populations morphology. The aim of the present study was to examine the 
morphological variation of A/burnus musselensis from all basins in Iran to evaluate the differences 
among its populations. 


MATERIAL AND METHODS 

In total 705 specimens of A. mossulensis were caught from 27 rivers of five basins in Iran, including 
Bushehr, Fars and Tigris basins (Karun, Karkheh and Diyala river basins) (Fig. 1, Table 1) in 2010 
by a seine net (5mm mesh sized), then the specimens were fixed in buffered 10% formalin after 
anesthetized in 1% clove oil solution at the sampling site and transferred to the laboratory for 
further examination. Twenty-two traditional morphometric characters were measured: total length, 
fork length, standard length, body height, head length, head width, snout length, eye diameter, cheek 
length, eyes distance, caudal peduncle length, caudal peduncle height, predorsal length, postdorsal 
length, preanal length, dorsal fin base length, dorsal fin height, anal fin base, length, anal fin height, 
ventral fin length, pectoral fin length and pectoral-ventral fin distance (Fig. 2). 

Distances were measured using digital caliper to the nearest 0.1 mm. As variation should be 
attributable to body shape differences and not related to the relative size of the fish; in order to 
eliminate any variation resulting from allometric growth, all morphometric measurements were 
standardized according to Elliott et al (1995): 

M4 = M @, / L)” 
Where M is the original morphometric measurement, M,,, the size adjusted measurement, L, the 
standard length of fish and L, is the overall mean of standard length for all fish from all samples for 
each variable. The parameter b was estimated for each character from the observed data as the slope 
of the regression of log M on log L,, using all specimens. The results derived from the allometric 
method were confirmed by testing significance of the correlation between transformed variables and 
standard length. Eleven meristic characters including pectoral fin soft ray, ventral fin soft ray, anal 
fin hard ray, anal fin soft ray, dorsal fin hard ray, dorsal fin soft ray, caudal fin soft ray, scales on 
lateral line, scales above lateral line, scales below lateral line, Circumpeduncle scales, predorsal scales, 
keel scales and gill arch spines, on the left side whenever possible were counted. Counts and 
measurements were done under a stereomicroscope. Also weight was measured with a Digital scale 
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to the nearest 0.1 gram. Univariate Analysis of Variance (ANOVA) was performed for each 
morphometric and meristic character to evaluate the significant difference between the five basins. 
Morphometric characters that had significant differences were used for principal component analysis 
(PCA), discriminant function analyses (DFA) and cluster analysis (CA). In the present study linear 
discriminant function analyses (DFA), principal component analysis (PCA) and cluster analysis (CA) 
were employed to discriminate the five populations. The Wilk's lambda was used to compare the 
difference among all groups. The DFA was used to calculate the percentage of correctly classified 
(PCC) fish. A cross-validation using PCC was done to estimate the expected actual error rates of the 
classification functions. As a complement to discriminant function analysis, morphometric and 
meristic distances between the populations of five basins were inferred to Cluster analysis (CA) 
(Veasey et a/., 2001). Statistical analyses were performed using the PAST version 2.17b, SPSS version 
19 software package and Excel 2013. 
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FIGURE 1. Iran map, basins and sampling sites (Based on Keivany ef a/., 2016d). 


76 IRANIAN JOURNAL OF ANIMAL BIOSYSTEMATICS Vol.14, No.2 


TABLE 1. Number, mean total length and weight of A. mossulensis from different basins of Iran. 


Basin River Counts 
Bushehr Fahlian 12 
Sheshpit 39 
Fars Safid 25 
Tang Boragh 13 
Zangiabad 40 
Karun Khersan 46 
Marbor 15 
Eivan Abbasi 30 
Bonestan 17 
Shirvan 11 
Dopolan 14 
Aligodarz 18 
Sepidar 15 
Darband 12 
Beshar 17 
Karkheh Gamasiab 14 
Sarab Gamasiab 20 
Malayer 62 
Kashkan 56 
Diyala Abdanan 11 
Haran 40 
Kangir 25 
Leileh 17 
Little Zab 16 
Rabat 86 
Shoshir 10 
Sirvan 24 
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FIGURE 2. A photo of Adburnus mossulensis with some morphometric measurements. 


MORPHOLOGICAL VARIATIONS OF ALBURNUS MOSSULENSIS 77 


RESULTS 

Morphometric characters 

The results of ANOVA for morphometric characteristics among five populations of A. mossulensis 
from the five basins of Iran are shown in Error! Reference source not found.. Significant 
differences (p<0.05) among populations of Bushehr, Fars, Karkheh, Karun and Diyala basins were 
observed for all the 22 morphometric characters (Error! Reference source not found.). These 
significant variables (except standard length and weight) were used for further multivariate analysis 
(PCA, DFA and CA). In order to determine which morphometric measurement most effectively 
differentiates populations, the contributions of variables to principal components (PC) were 
examined. Being KMO coefficients approximately more than 0.6 indicate that PCA method will 
suitable for our data (Kaiser, 1974). For morphometric characters, the KMO coefficient were 
obtained as 0.946 and for morphometric characters that is explaining of appropriation of this test at 
good and medial level, and the Bartlett’s Test of sphericity is significant (PS0.01) (Error! Reference 
source not found.). In PCA analysis, the characters with an eigenvalues more than 1 were included 
and others discarded (Error! Reference source not found.). 

The Wilks’ lambda test of discriminant analysis indicated significant differences in 
morphometric characters of the five populations from basins. In this test, all functions were highly 
significant (PS0.01) ( 

). Nimalathasan (2009) mentioned out that factor loading greater than 0.30 are considered 
significant, 0.40 are considered more important and 0.50 or greater are considered very significant. 
Principal component analysis of 22 morphometric characters showed that PC 1 accounts for 
68.556% of the variation, PC II for 13.274% and PC HI for 7.395% (Error! Reference source not 
found.) and that the most significant weightings on PC I were from B.H, H.L, H.W, Sn.L, E.D, 
Po.L, E.W, C.L, C.H, Pr.D.L, Po.D.L, D.H, A.B, A.H, V.L, P.L and P-V.L, on PC II were from T.L, 
P.L, B.H, H.W, E.W, C.H and on PC III were from Pr.A.L and D.B (Error! Reference source not 
found.). The rotated (Varimax) component loadings for the three components (factors) are 
presented in Error! Reference source not found.. 

For the discriminant analysis, the averages of percentage of correctly classified (PCC) were 
80.5% for morphometric characters. High classification success rates were obtained for Bushehr 
(72.5%), Fars (98.7%), Karkheh (66.7%), Karun (100%) and Diyala (64.6%) stocks indicating a high 
correct classification of individuals into their original populations with respect to morphometric 
characters (Error! Reference source not found.). Error! Reference source not found. indicates 
the coordinates of the five populations in the two first axes of DFA. In this analysis there was a low 
degree of separation among three populations from Bushehr, Karkheh and Diyala basins and high 
degree of separation among Fars and Karun basins population with selves and with other 
populations of A. mossulensis. 

Cluster analysis based on Euclidean distances in morphometric characters among the groups 
of centroids in A. musselensis populations resulted into grouping of populations in three clusters, the 
Fars basin, the Karun basin and the other remaining three basins (Error! Reference source not 
found.). This grouping pattern was similar to discriminant functions analysis (DFA) in Error! 
Reference source not found.. 


Meristic characters 

The results of ANOVA for meristic characteristics among five populations of A. mossulensis from the 
five basins of Iran are shown in Error! Reference source not found.. Significant differences 
(p<0.05) among populations of Bushehr, Fars, Karkheh, Karun and Diyala basins were observed for 
all the meristic characters except counts of anal and dorsal fin spine (Error! Reference source not 
found.). So, these significant variables were used for further multivariate analysis (PCA, DFA and 
CA). In order to determine which meristic measurements most effectively differentiate populations, 
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the contributions of variables to principal components (PC) were examined. Being KMO 
coefficients approximately more than 0.6 indicate that PCA method will be suitable for our data 
(Kaiser, 1974). For meristic characters, the KMO coefficient were obtained as 0.722 and for meristic 
characters which explains the appropriation of this test at good and medial level, and the Bartlett’s 
Test of sphericity was significant (PS0.01) (Error! Reference source not found.). 

TABLE 2. The mean+SD and results of ANOVA for morphometric characters of A. mossulensis from the five 
basins of Iran. 


N Bushehr Fars Karkheh Karun Diyala 

o. Abb. Characters 51 78 152 194 229 

1 TL Total Length 89,.1+13.2b  77.3426,74 82.0+20.4¢ — 85.3424.1¢ — 73.9427.44 
2 FL Fork Length 80.8412.2> —70,.2+24.6¢ 75.0+19.1¢  79.0+20.9* — 67.5425,44 
3 SL Standard Length 73.5+11.4  63.3422.9 67.4417.4  71.7£17.4 — 60.9423.5 
4 BAH Body Height 16.642.74 = 13.15.38 14.24+3,9> 16.445,3¢ 12.8+5.53 
5 HL Head Length 18.7+£2.8¢ 17.0£5.6¢ 16.8+4.0° 17.644.54 15.245.4+ 
6 H.W Head Width 8.9+1.4¢ 7.5£2.6% 7.7£1,9% 8.5£2.3> 7.2£2.68 

7 Sn.L Snout Length 5.0£0.8¢ 4.441,2> 4,241.1» 4.441.3> 3,841.58 

8 ED Eye Diameter 5.1£0.6¢ 4.34115 4.6£0.9> 4441.08 4.041.3+ 

9 PoL Cheek Length 9.641.5> 8.83.25 8.5+2.28 9,242.48 1942.99 
10 EW Eyes Distance 5.6£0.8¢ 4,541.68 4.941.2> 5.241.4> 4.441.8* 
11 CL Caudal Length 18.542.9 15.65.48 16.844,3> 17.845.08 15.245.8* 
12 CH Caudal Height 7.341.1¢ 5.7£2.48 6.3+1.7> 7142.16 5.7£2.4ab 
13, Pr.D.L  Predorsal length 38.1+6.0¢ — 33,7£12.1¢4 35.0+£9.3¢ 37.8411.1% — 31.3+11.9> 
14 Po.D.L  PostDorsal Length 29.044.5¢ 23.88.98 25.9+6.8> 27.2+8.58 23,549.38 
15 PrA.L  PreAnal Length 47,847.64  41.6415.48 44.1411.54  47.3414.1> — 39.8415,.2¢ 
16 DB Dorsal fin Base length 8.1+1.4> 6.842.3¢ 7.3£2.0% 7,742.48 6.7+2.54 
17. D.H Dorsal fin Height 15.042.24 = 12.5+3.8> 13.543.2¢ 13:143.34 12.144.2> 
18 AB Anal fin Base length 10.0£1.6° 8143.14 9.0+2.4> 8.842.68 8.0+3.08 
19 AA Anal fin Height 11.841.7° — 10.0+3.0¢ 10.8£2.5'¢ — 10.5+2.64 9,5+3,2> 
20. V.L Ventral fin Length 11.441.8° 9,643.63 10.4+2.7> 10,242.94 913.5% 
21 P.L Pectoral fin length 14.842,3¢ 11.944.4> 13.143.4% — 13,0+3,5> 11.6+4.28 
22 P-V.L Pactoral-Ventral Distance 16.942.7> = 13,945.48 15.544.4> 17.6£5.8¢ 14.145.6> 
2 W Weight 6.443,3> 4.845.4* 4,743.48 TALT.8P 4,744.74 


No.: Numbers, Abb.: Abbreviation. a, b, c, d and e are result of Duncan grouping. 


TABLE 3. Eigenvalues, percentage of variance, percentage of cumulative variance and Canonical Correlation 
for the Canonical Discriminant Functions in case of morphometric characters of A. mossulensis from the five 
basins of Iran. 


Function Eigenvalue % of Variance Cumulative % Canonical Correlation 
1 65.964 82.5 82.5 0.993 
2 13.456 16.8 99.3 0.965 
3 0.378 0.5 99.8 0.524 
4 0.183 0.2 100.0 0.393 


TABLE 4. Result of Wilks' lambda test for verifying difference among five populations of A. mossulensis from 
the five basins of Iran. When morphological measurements are separately compared using discriminant 
function analysis. 


Test of Function(s) Wilks’ Lambda Chi-square Df Sig. 

1 through 4 0.001 5081.079 84 0.000 
2 through 4 0.042 2180.215 60 0.000 
3 through 4 0.613 337.143 38 0.000 


4 0.845 115.877 18 0.000 
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FIGURE 3. Coordinate plot of A. mossulensis from the five basins of Iran according to the first two 
discriminant functions from morphometric data analysis. 


TABLE 5. KMO test and Bartlett's Test in case of morphometric characters for of A. mossulensis from the five 


basins of Iran. 


Kaiser-Meyer-Olkin (KMO) Measure of Sampling Adequacy. 0.946 
Bartlett's Test of Sphericity Approx. Chi-Square 28692.264 
Df 210 
Sig. 0.000 


TABLE 6. Eigenvalues, percentage of variance and percentage of cumulative variance for the principal 


components in case of morphometric characters for of A. mossulensis from the five basins of Iran. 


Component Eigenvalues % of Variance Cumulative % 
1 14.397 68.556 68.556 
2 2.788 13.274 81.830 
3 1.553 7.395 89.226 


TABLE 7. Factor loadings for the principal components and correlations between the morphometric 
characters for A. mossulensis from the five basins of Iran. 


PC 1 PC 2 PC 3 PC 1 PC 2 PC 3 
Adj-T.L -0.931 Adj.Pr.D.L 0.904 
Adj-F.L -0.911 Adj.Po.D.L 0.898 
Adj.B.H 0.818 0.498 Adj.Pr.A.L -0.949 
Adj.H.L 0.897 Adj.D.B 0.959 
Adj.H.W 0.871 0.405 Adj.D.H 0.919 
Adj.Sn.L 0.891 Adj.A.B 0.882 
Adj.E.D 0.886 Adj.H 0.909 
Adj.Po.L 0.888 Adj.V.L 0.959 
Adj.E.W 0.877 0.421 Adj.P.L 0.949 
Adj.C.L 0.819 Adj.P-V.1 0.794 0.451 
Adj.C.H 0.845 0.492 
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Adj: adjusted 


TABLE 8. Percentage of specimens classified in each group and after cross validation for 
morphometric characters for A. mossulensis from the five basins of Iran. 


basins Bushehr Fars Karkheh Karun Diyala 
Original % 

Bushehr 76.5 0.0 19.6 3.9 
Fars 98.7 13 
Karkheh 15.8 0.0 69.1 0.7 14.5 
Karun 100.0 

Diyala 13:5 20.1 66.4 
Cross-validated % 

Bushehr 72.5 21.6 5.9 
Fars 98.7 13 
Karkheh 16.4 66.4 0.7 16.4 
Karun 100.0 

Diyala 14.8 20.5 64.6 


Fars 
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FIGURE 4. Dendrogram derived from cluster analyses of 21 morphometric characters on the basis 
of Euclidean distance for A. musselensis populations. 


The Wilks’ lambda tests of discriminant analysis indicated significant differences in metistic 
characters of the five populations. In this test, all functions were highly significant (PS0.01) (Error! 
Reference source not found.). In PCA analysis the characters with an eigenvalues more than 1 
were included and others discarded (Error! Reference source not found.). Principal component 
analysis of 14 meristic characters showed that PC1 accounts for 21.23% of the variation, PC II for 
17.87%, PCIII for 9.39% and PC IV for 8.82% (Error! Reference source not found.3) and that 
the most significant weightings on PC I were from scales on lateral line, scales above lateral line, 
Circamucaudal scales and Predorsal scales, on PCII were from anal fin soft ray, pectoral fin soft ray, 
ventral fin soft ray and keel scales, on PC III were from caudal fin soft ray and gill arch spines and 
on PCIV were from dorsal fin soft ray, scales down lateral line and keel scales (Error! Reference 
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source not found.). The rotated (Varimax) component loadings for the three components (factors) 
are presented in Error! Reference source not found.. 

For the Discriminant Factor Analysis, the averages of percentage of correctly classified 
(PCC) were 93.0% original data and 92.3% for Cross-validated data for meristic characters. High 
classification success rates were obtained for Bushehr (100%), Fars (100%), Karkheh (92.9%), 
Karun (86.2%) and Diyala (95.1%) stocks indicating a high correct classification of individuals into 
their original populations with respect to meristic characters (Error! Reference source not found.). 
Error! Reference source not found. indicates the coordinates of the five populations in the two 
first axes of DFA. In this analysis there was a low degree of separation among three populations 
from Bushehr, Karkheh and Diyala basins and high degree of separation among Fars and Karun 
basins population with selves and with other populations of A. mossulensis. Clustering analysis based 
on Euclidean distances in metistic characters among the groups of centroids in A. musselensis 
populations resulted into grouping of populations, first, Bushehr basin separated, then Karkheh and 
other three basins were in one group (Fig. 6). 


TABLE 9. The results of ANOVA for morphometric characters of A. mossulensis from the five basin 
of Iran. 


Pectoral fin soft ray Ventral fin soft ray 
0 % Frequency of each 
% Frequency of each count 
count 

basins meant+SD 13. 14 15 16 17:18 mean+SD 7 8 9 10 
Bushehr 15.2+0.8> 0 18 47 29 6 0 8.6£0.5> 2 37 61 O 
Fars 15.6+1.1¢ 0 15 3331 17. 4 8.7£0.5> 0 33 67 O 
Karun 14.6£0.8¢ 21 34 15 11 3 1 8.340.22 2 70 28 O 
Karkheh 15.3+1.1b¢ 0 3.3, 23 34.9 25 13.8  8.7£0.6> 0 29 71 +O 
Diyala 15.0+1.5> 14 28 = =28 ~— 13 7 10 8.6+0.9> 2 41 #56 1 


TABLE 9 Continued. The results of ANOVA for morphometric characters of A. mossulensis from 
the five basin of Iran. 


Anal fin hard ray Anal fin soft ray 

% Frequency of % Frequency of each count 

each count 
basins meantSD 2 3 meantSD 9 10 11 12 13 
Bushehr 3.00.0? 0 100 11.3241.5¢ 0 5.9 41.2 47 5.9 
Fars 3.00.02 0 100 10.9+0.8> 3 28 49 19 1 
Karun 3.00.08 1 99 14.60.82 21 36 16 11 1 
Karkhe 3,040.12 0.7 99.3 11.140.6>¢ 0 11.8 64.5 21.7 2.0 
Diyala 3.00.08 0 100 11.040.7» 1 26.5 23.4 37 7 


TABLE 9 Continued. The results of ANOVA for morphometric characters of A. mossulensis from 
the five basin of Iran. 


Dorsal fin hard ray Dorsal fin soft ray 
0 z 0 ; 

nen oD % Frequency of each count anaes % Frequency of each count 
basins 2 3 7 8 9 
Bushehr 2.7£0.54 15 85 7.940.4a> 2, 14 84 
Fars 3.00.08 0 100 7.90.3» 8 1 91 
Karun 3,040.1? 3 97 7.940.125 12 85 S) 
Karkhe 3.00.02 4.6 95.4 8.040,.2> 4 95 1 


Diyala 3.00.12 1 99 7.80.42 21 78 1 
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TABLE 9 Continued. The results of ANOVA for morphometric characters of A. mossulensis from 


the five basin of Iran. 


Scal pears Circamu 
Scales on ey down : Predorsal Keel Gill arch 
Caudal fin soft ray : lateral caudal : 
lateral line } lateral scales scales spines 
line . scales 
line 
Frequency of 
each count (%) 
basins meantsd 18 19 20 meantsd meantsd meantsd meantsd meantsd meantsd meantsd 
Bushehr 18.90.42 10 86 4  64.2+3.48 12.741.0* 4.940.427 19.841.98  30.242.5° 9,.841.6¢ 11.30.84 
Fars 19,040.35 3 94 4  71.244,3¢ 14.740.8¢ 5340.8 21.842.95 35.142.1¢ 10.041.4¢ 13.441.3¢ 
Karun 18.80.12 15 85 O 65.7£13.85 13.342.5> 4.940.9° 20.443.9> 31.146.1>5 7.541.484 11.7+0.5> 
Karkheh 19.040.2 2.6 954 2 75.044.84 14.840.6¢ 5.2+0.5>* 22,041.2> 34.442.8¢ 9,241.25 12.14£0.9¢ 
Diyala 18.940.44 11 88 1 69.9+4.9b¢ 14,340.95 5.140.5> 21.341.0> 33.643.1> 8.941.1> 12.641.14 


TABLE 10. Eigenvalues, percentage of variance, percentage of cumulative variance and Canonical 
Correlation for the Canonical Discriminant Functions in case of morphometric variables of A. 
mossulensis from five basins of Iran. 


Function Eigenvalue % of Variance Cumulative % Canonical Correlation 
1 18.2314 84.3 84.3 0.974 
2 2.1028 9.7 94.0 0.823 
3 0.9532 4.4 98.4 0.699 
4 0.3428 1.6 100.0 0.505 


TABLE 11. Result of Wilks' lambda test for verifying difference among five populations of A. 
mossulensis from the five basins of Iran. When morphological measurements ate separately compared 
using discriminant Function analysis. 


Test of Function(s) 
1 through 4 

2 through 4 

3 through 4 

4 


Wilks' Lambda 


0.006 
0.123 
0.381 
0.745 


Chi-square 


2895.057 
1200.955 


552.27 


1 


168.586 


df 
68 
48 
30 
14 


Sig. 

0.000 
0.000 
0.000 
0.000 


TABLE 12. KMO test and Bartlett's Test in case of morphometric variables for of A. mossulensis from 


the five basin of Iran. 


Kaiser-Meyer-Olkin (KMO) Measure of Sampling Adequacy. 0.722 
Bartlett's Test of Sphericity Approx. Chi-Square 1161.644 
df 66 
Sig. 0.000 
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FIGURE 5. Coordinate plot of A. mossulensis from the five basins of Iran according to the first two 
discriminant functions from morphometric data analysis (1: Bushehr, 2: Fars, 3: Karkheh, 4: Karun, 


5: Diyala). 


TABLE 13. Eigenvalues, percentage of variance and percentage of cumulative variance for the 
principal components in case of morphometric variables for of A. mossulensis from the five basins of 


Tran. 
Component Eigenvalues % of Variance Cumulative % 
1 2.547 21.228 21.228 
2 2.145 17.873 39.101 
3 1.127 9.392 48.493 
4 1.058 8.817 57.310 


TABLE 14. Factor loadings for the principal components and correlations between the 
morphometric variables for A. mossulensis from the five basins of Iran. 


Component Component 

1 2 3 4 1 2 3 4 
Dorsal fin soft ray -0.685 Scales up lateral line 0.739 
Anal fin soft ray 0.538 Scales down lateral line 0.598 
Pectoral fin soft ray 0.709 Circamucaudal scales (0.794 
Ventral fin soft ray 0.697 Predorsal scales 0.709 
Caudal fin soft ray 0.775 Keel scales 0.701 0.415 
Scales on lateral line 0.752 Gill arch spines 0.713 
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TABLE 15. Percentage of specimens classified in each group and after cross validation for 
morphometric characters for A. mossulensis from the five basins of Iran. 


Bushehr Fars Karun Karkheh Diyala 
Original (93.0%) 
Bushehr 100.0 
Fars 100.0 
Karkheh 93.6 6.4 
Karun 87.0 13.0 
Diyala 0.6 4.3 95.1 


Cross-validated (92.3%) 
Bushehr 100.0 


Fars 100.0 

Karkheh 92.9 7A 
Karun 86.2 13.8 
Diyala 4.9 95.1 


Bushehr 
Karkheh 


Fars 


Diyala 


Karun 


FIGURE 6. Dendrogram derived from cluster analyses of 12 meristic characters on the basis of 
Euclidean distance for A. musselensis populations. 


DISSCUSSION 

The multivariate analysis of morphometric and meristic characters in this work indicates the 
existence of morphologically differentiated groups of A. mossulensis and classified populations along 
the basins into five distinct groups. Body size is a major source of error in statistical analyses of data 
and plays a major role in morphometric analysis ([zeng, 2004); however, in this study, allometric 
transformation was successfully removed the body size effect, so all significant variations showed 
the body shape differences when it was tested using ANOVA and other multivariate analysis. 

In the present study, highly significant morphological variations were detected among A. 
mossulensis populations from Bushehr, Fars, Karkheh, Karun and Diyala. The results obtained from 
ANOVA analysis in our study showed that all of 21 transformed morphometric and 12 meristic 
characters were significantly different in populations of A. mossulensis existing in five basins of Iran 
that exhibits a high phenotypic variation among these populations. 

Discriminant Function Analysis (DFA) is a useful method to distinguish different stocks of a 
same species, concern to stock management programs (Karakousis ef a/, 1991). The results of DFA 
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confirmed a high differentiation among the five populations of A. mossulensis in the study areas. This 
separation pattern was greatly confirmed by PCA, where showed that the populations were distinct 
from each other. 

The reasons of morphological differences between populations are often quite difficult to 
explain (Bookstein, 1991; Poulet e¢ a/., 2004). It has been suggested that the morphological 
characteristics of fish are determined by genetic, environment and the interaction between them 
(Swain & Foote, 1999; Pinheiro ef a/, 2005; Anvarifar et a/, 2010; Keivany ef a/, 2016b). The 
population from Fars basin (in the central part of the Iran) and Karun basin (in the southwest part 
of the Iran) are more distinct from each other and other populations, on the other, populations from 
Bushehr, Karkheh and Diyala basins were overlapping. It is well known that morphological 
characteristics can show high plasticity in response to differences in envitonmental conditions. This 
raises the probability that phenotypic plasticity may itself be adaptive, allowing fish populations to 
change their appearance to match their ecological environments (Swain & Foote, 1999; Mousavi, 
2011; Zamani-Faradonbe ef a/, 2015b). The environmental factors that are dominant during the early 
development stages, when individual’s phenotype is more amenable to environmental influence, is 
of particular importance (Pinheiro e¢ a/., 2005). 

Therefore, the distinctive environmental conditions in divers rivers that flowing in Bushehr, 
Fars, Karkheh, Karun and Diyala basins, might cause the morphological differentiation between 
these five populations. The phenotypic variability may not necessarily reflect population 
differentiation at molecular level (Ihssen e¢ a/., 1981). In overall, fishes reveal greater difference in 
morphological characters both within and between populations than other vertebrates, and are more 
susceptible morphological variations that were effected by environmental conditions (Thompson, 
1991; Wimberger, 1992; Turan ef a/, 2006; Zamani-Faradonbe ef a/., 2015c), so these difference 
might reveal different physicochemical conditions, diets, prey types, food availabilities, feeding 
pattern or other features (Rezaei ef a/., 2012). 

The present study proposes high morphological differentiation among A. mossulensis 
populations along the five basins in Iran. The results also suggest these morphologically different 
populations should be considered as distinct stocks in the basins in fisheries management and 
commercial exploitation of this species and any stock enhancement program. Nevertheless, future 
studies on determination of population structure will be elucidated using biochemical and molecular 
genetics methods. 
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